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Stability of solitary waves in a thin inextensible and unshearable rod of infinite 
length is studied. Solitary-wave profile of the elastica of such a rod without torsion 
has the form of a planar loop and its speed depends on a tension in the rod. The 
linear instability of a solitary-wave profile subject to perturbations escaping from 
the plane of the loop is established for a certain range of solitary-wave speeds. It is 
done using the properties of the Evans function, an analytic function on the right 
complex half-plane, that has zeroes if and only if there exist the unstable modes of 
the linearization around a solitary-wave solution. The result follows from compar- 
ison of the behaviour of the Evans function in some neighbourhood of the origin 
with its asymptotic at infinity. The explicit computation of the leading coefficient 
of the Taylor series of the Evans function near the origin is performed by means of 
the symbolic computer language. 
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1. Introduction 

Dynamics of inextensible and unshearable thin rods of the infinite length described 
in the framework of the Bernoulli-Euler beam model is governed by the following 
equations (in dimensionless form) 

t« = (*»-')« + ie - 

TV; = 1, 

n-»i, t- 2i3 ->o, £->±oo, (1.1) 

where £ is the dimensionless arc- length along the elastica, = t 1 , i = 1,2,3 are 
the components of the unit vector tangent to the elastica and p = P — 1, P being 
the absolute value of tension in the rod; summation is assumed under repeating 
indices and subscripts t and £ denote differentiation with respect to these variables. 
Hereafter, the different notations for the same components are given both with lower 
as well as upper indices; it is done for the sake of convenience (when a component 
already has the subscripts denoting differentiation or some other) and also in the 
case when the rule about summation under repeating indices is assumed. 

The classification of the forms of elastica was for the first time presented by 
Euler, who derived the ordinary differential equation of elastica (Love [19]). For the 
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Figure 1. Solitary- wave profile on an elastica and tangent vector {t° , T2,0} t . 

elastica of the infinite length its profile admits the configuration described by the 
solitary-wave solution of (1.1), having the form of the planar loop (fig.l): 



p = -p° = -6(1 - c 2 ) sechVl ~c 2 C, t 3 = 0, n = if = 1 - 2scchVl - c 2 C, 



T2 



= r 2 ° = -2sech 2 y/\ - c 2 Csinh\/l - c 2 C, £ = £ - ct, < c < 1. 



(1.2) 



The formula (1.2) gives the particular localized solution, where the tension P is 
distributed in the rod in a following way: it increases from the minimal value at 
the top of the loop (for 1 — c 2 > 1/6, for example, this value is negative, that 
corresponds to the compression of the rod at the place of the loop's localization) 
to the maximal value at both infinities. 

The system of equations (1.1) is a hamiltonian one and it can be written as 



i _ d_SH_ l _^_$Ji 



l. 



where 



H= j [( V y+T tl :Tl+P(nT*-l)]dC 

— oo 

Along with the hamiltonian H one has two more conserved quantities: 



Q 



oo oo 

yv-r^de, i?=yv-i)^, ri=\, t^=o. 



In this paper we prove linear instability of certain solitary waves with respect 
to transverse perturbations <5t 3 ^ 0. It is shown that there exists exponentially 
growing with time eigenfunction being the solution of the linearization about a 
solitary wave when its dimensionless speed does not exceed 1/^/2. The existence 
of the unstable eigenfunction is expected to imply the physical instability, i. e. the 
decay of the loop under the transverse perturbations. The opposite result about 
the stability, when the loop is not destroyed, takes place for plane perturbations 
5t 3 = 0, Sti, 2 7^ 0. Mathematically this result is expressed by the Theorem 1.1 
below. 
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Let us denote (j) c = {rf - 1, t 2 °, 0, v%, 0} T (uj = -c^ - 1), u§ = -ct 2 °), 
w(t) = {ti - l,7vj,7^,«i,«2,«3} T and let X = H^R) x fT^K) x H^R) x i 2 (M) x 
L 2 (M) x i 2 (R) be a Hilbert space, || • || is the norm in X. In Beliaev and Il'ichev [6] 
using the ideas from Grillakis et al. [12] the following theorem is proved relating to 
the nonlinear stability of solitary waves (1.2) subject to perturbations of rf 2 only. 

Theorem 1.1. Let t 3 = and for any given Wq e X near <p c in X , ||wo — <f> c \\ < 1, 
there exist T = T("f) > and a vector function w(t) G C([0,T), X) (continuous 
with values in X) such that w(0) = w = {tq, 0, v l , 0} T , I = 1,2, and for all 
0<t<T, r^roi = r t T l = 1 and H(w) = H(w ), Q(w) = Q(w ), J2(w) = J2(w ). 
T/ien, /or all e > there exists S > swc/i i/iai i/ ||w — C || < (5, i/ien 

sup inf ||w(t) - ^> c (- + w)|| < e. 
t>0 weR 

The theorem was proved by demonstrating that the set of translates of a solitary 
wave gives a local minimum to the invariant functional 

F(w) = H(w) + cQ(w) - i?(w) 

on the closed submanifold Tj-r* = 1 in X. This, in its turn, follows from the fact that 
the 'linearized hamiltonian' H(<j) c ) = S 2 F{<p c )/5 2 w has exactly one zero eigenvalue 
and positive spectrum bounded away from zero. The analysis of Beliaev and Il'ichev 
[6] is extended in Dichmann et al. [9] to the case when the rotational kinetic energy 
(which is small compared to the total energy for thin rods) is not neglected and 
the same stability result for subsonic solitary waves was established. For St^ — 
= t 7^ 0, however, the operator Ti{<j> c ) has an extra negative eigenvalue and the 
sufficient conditions for the stability of bound states as they appear in Grillakis et 
al. [12] are violated. The system of equations (1.1) is also invariant under group of 
rotations (see sec. 6). However, the stability theory given in Grillakis et al. [13] for 
hamiltonian systems with several symmetries is not applicable in the case under 
consideration by the reason that solitary wave orbit is in fact the set of translates of 
a solitary wave, and therefore the conserved quantity associated with the rotational 
symmetry, being nonlocal, identically equals zero at a solitary wave. 

The linearized equations (1.1) for 0Y1.2 and r decouple and we analyze here only 
the equation for r since the stability result in the planar case is known. A solitary 
wave is called linearly unstable if there exist solutions of the linearized equation of 
the form 

T = r(t,t) =e xt w(\X) (1.3) 

with w(X, () decaying exponentially as £ — > ±00 and Re A > 0. The function w(X, () 
obeys the eigenvalue problem 

The equation (1.4) has the same form as the corresponding equation for the stability 
problem discussed in Alexander and Sachs [2]. That is why, the theoretical results 
of this paper are also valid for our analysis. Yet, the coefficient p° in (1.4) does 
not coincide as a function of £ with the corresponding coefficient in Alexander and 
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Sachs [2] and, as a consequence, our computations and the result are different from 
the computations and the result of these authors. 

We get the instability result by means of the Evans function. The Evans function 
D(X) is constructed as an analytic on the right complex half-plane function of the 
spectral parameter; D(X) has zeroes if and only if there exist the unstable modes 
of the linearization around a solitary-wave solution. The use of an Evans function 
to get instability result for solitary waves appears in Pego and Weinstein [20]. In 
parabolic problems, the ideas of Evans [10] were further developed by Jones [15] 
and Alexander et al. [3] (see also Kapitula [16] and references therein). Alexander 
and Sachs [2] extended the technique of Pego and Weinstein [20] to the case of the 
Bousssinesq-type equations where two decaying modes are present in the model. 
The application of the Evans function method for instability problems in various 
fields of hydrodynamics and physics can be found in Swinton and Eglin [22], Pego 
et al. [21], Alexander et al. [4], Gardner and Zumburin [11], Kapitula [17], Kapitula 
and Sandstede [18], Afendikov and Bridges [1], Bridges et al. [8]. 

The Evans function is real on the real axis and tends to unity as |A| — > oo. Our 
instability result follows from comparison of the behaviour of the Evans function 
in some neighbourhood of the origin with its asymptotic at infinity. We show that 
in a small neighbourhood of the origin one has D(X) < for real positive A. More 
precisely, the following theorem holds. 

Theorem 1.2. The Evans function D(X) constructed for the eigenfunction w(X, () 
from (1.3) is analytic in a neighbourhood of the origin and its converging Taylor 
series in this neighbourhood reads 

1 _ o r 2 °° 
^) = - 4(1 _ c2)3 A 2 + ^e„(c)A". (1.5) 

^ ' n=3 

It follows then, that for c 2 < 1/2 the function D(X) has a zero somewhere at 
A = Ao on the positive real axis, that, in its turn, implies the existence of the 
unstable eigenfunction w(Xo,() from (1.3). 

In section 2 we give the formulation and linear stability problem. In section 3 
we discuss analytic properties of the solutions of the linearized system as A varies 
in a zero neighbourhood. In section 4 the Evans function is introduced and its 
behaviour for large |A| is described. In section 5 we describe computations (made 
by means of the symbolic language Mathematica 4.0) of the leading coefficient of 
the Taylor series of D(X) in a neigbourhood of the origin. The explicit form of key 
expressions in our computations are listed in Appendix A. In section 6 we discuss 
our conclusions. 

2. Formulation and linear stability problem 

We consider a thin inextensible and unshearablc elastic rod of infinite length initially 
coinciding with the X\ axis of a Cartesian coordinate system. The total energy of 
such a rod is the sum of the kinetic energy and the bending energy, the torsion 
energy is neglected. The respective energy densities K and W are given by the 
expressions 

K = -pSx\x iu W = -I Excite , 
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where Xi, i — 1, 2, 3 are the coordinates of the points on a neutral line (elastica) of 
the rod, p is the density of the rod, S is the area of the cross section of the rod, pi 
is the moment of inertia of the cross section of the rod about the line, orthogonal 
to the principal plane of bending X\x 2 , E is the Young module, £ is the arc-length 
along the elastica. The elastica is given by the equations x l = x l (t;, t). For thin rods 
the rotational part of the kinetic energy is small in comparison with the kinetic 
energy of points of the elastica [6] and we neglect it here. 

The equations of motion can be obtained by taking the variations of the la- 
grangian A, 



A = I 
2 



(pSx\xit — IEx££Xi££)d£ dt 



to 



under the condition of incxtensibility 



Xi£ x^ 



1. 



These equations read 



where P(£, t) = p(£, t) +Poc is the Lagrange multiplier, p 
the scaling transformations in (2.1) 



as £ 



(2.1) 
±oo. Making 



P -» PooP, £ -» y/lE/poot, t -> ^pSIE/p^t, 



x k . 



and preserving the old notations one gets (1.1) with T{ 

Linearizing the equations (1.1) about the solitary-wave solution (1.2) we get 
the following equation for the perturbation St 3 = r of the third component of the 
tangent vector (1.2) 



Ttt = -(P°r)«+r« -t««. 



(2.2) 



We seek the solution of (2.2) in the form (1.3). After substitution of (1.3) into the 
equation (2.2) it transforms into the equation (1.4). The last equation, in its turn, 
can be written in the matrix form 



y' = M(A,C)y, 

y = {yi, 2/2,2/3, 2/4> T , 2/i=io, y2=w', y 3 =w", y 4 = w" 
where prime denotes differentiation with respect to (, and 



(2.3) 



M(X,C) 
The adjoint equation reads 



V 



1 \ 
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A 2 - p°" 2Ac - 2p ' l-c 2 -p° J 



(2.4) 
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or in the matrix form 

z' = -zM(\,(), (2.5) 

z = {zi,z 2 , z 3 , z 4 }, z 4 = to*, z 3 = -w* , z 2 = w* - (1 - c 2 - p°)w*, 
z x = -w*'" + (l-c 2 -p°)w*' - (2Xc-p°')w*. 

3. Analytic properties of the solutions of the linearized 
system for A in the zero neighbourhood 

In this section we enumerate the results of Alexander and Sachs [2] concerning the 
analytic properties of the solutions of systems having the form of (2.3), (2.5) when 
A locates close to zero. 

Let us denote .Moo (A) = lim A4(\,(), M<*(A) (a = 1,2,3,4) the eigenvalues 

C— >±oo 

of _Moo(A), r Q (A) the right and l a (A) the left eigenvectors of the matrix .Moo(A). 
The matrix .Moo (A) is given by the same expression as the corresponding matrix in 
Alexander and Sachs [2] . For A = the vectors labeled by the subscript 4 represent 
the generalized eigenvectors of .Moo (A) (see formulas (3.9), (3.10) in Alexander and 
Sachs [2]). The eigenvectors are normalized so, that the first (for right eigenvectors) 
and last (for left eigenvectors) components are unity. 

The characteristic equation det [.Moo (A) — = 0, where £ denotes the unit 
matrix, has the form 

/i 4 - (1 - (?)[? — 2fic\ + A 2 = 0. (3.1) 
In Alexander and Sachs [2] the following lemma is proved. 

Lemma 3.1. For Re A ^ the equation (3.1) has two roots in each complex half- 
plane. 

Following [2] we denote Hi(X), /U3(A) the roots, lying in the left complex half- 
plane for Re A > (| Re/ii(A)| > | Re/i3(A)| for A in a zero neighbourhood). In the 
neighborhood of the origin the roots /Ufc(A), k = 1, 3 have the asymptotic form [2] 

Al3 (A) = --^ + 0(A 3 ). (3.2) 
1 — c 

The eigenvalues /ii.3(A) and the eigenvectors ri^A), L^A) may lose the analyticity 
property in the points where the eigenvalues are not simple any more or, in other 
words, when the corresponding roots of equation (3.1) are multiple. The roots of 
(3.1) become multiple when the resultant of the polynomial in the left hand side of 
(3.1) and its derivative equals zero. The resultant is 



16A 2 [(1 - c 2 ) 3 + (-8 - 20c 2 + c 4 )A 2 + 16A 4 ], 
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Figure 2. Curves A versus c where the eigenvalues are multiple. Along the branch A = A3 
the eigenvalues ^2 and /14 coincide; along the branch A = A5 the eigenvalues hi and /13 
coincide; along the branch A = A 6 the eigenvalues 1x3 and fj,2 coincide and along the branch 
A = A4 the eigenvalues /xi and /14 coincide. 

and it equals zero when A 1.2 = 0, and 



A 3 , 4 = ± 

^5,6 = ± 



a/8 + 20c 2 


-c 4 + 


c(8 + c 2 ) 3 / 2 




\/8 + 20c 2 


-c 4 - 


c(8 + c 2 ) 3 / 2 


4\/2 



(3.3) 



The nonzero branches (3.3) are pictured in Fig. 2. It is seen that for nonzero A the 
points where the eigenvalues are multiple are separated from zero when < c < 1. 
It is also seen that when A crosses the imaginary axis the eigenvalues /z 3 and /J4, 
coinciding with zero at A = 0, escaping from their complex half-planes. 

At A = one has /Lt 3 (0) = /^(O) = 0. Yet, the kernel of -Moo(A) — ^ 3 (A)£ 
spanned by r 3 (A) (13(A)) is uniformly one dimensional and the vectors themselves 
are continuous at A = (see formulas (3.9), (3.10) in Alexander and Sachs [2]). 
Hence, in a small enough neighbourhood of the origin, not including the points 
given by (3.3), the normalized vectors i"i ;3 (A) and li, 3 (A) are analytic. 

It follows from general theory of the ordinary differential equations that there 
exist solutions of (2.3), (2.5), satisfying 

lim c-^ A V fc (A,C) =r fc (A), 

lim e^z k {X,Q =lfe(A), fc = l,3. (3.4) 

00 

Analyticity of i"i >3 (A) li, 3 (A) implies the following theorem, proved in Pego and 
Weinstein [20] and generalized for the two mode case in Alexander and Sachs [2]. 

Theorem 3.2. Solutions yi, 3 (A, (), z i, 3 (A, C) are analytic in some neighborhood of 
A = 0. 
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4. Evans function 

Consider the vector fields y A (A, Q, z A (A, Q with the components 

Va/\t3 = Vi*yzf3 - yipysa, z£ A/} = zi a z 3 p - zi z 3a , a>0, a, = 1,2,3,4, (4.1) 

where yk a , Zk a are the components of the vectors yfc, Zk, correspondingly. We 
number the pairs a A [3 by the following way: 1 A 2 ^ 1, 1 A 3 ^ 2, 1 A4 ^ 3, 
2A3^4, 2A4^5, 3A4^6. The vectors y A (A, Q, z A (A, C) obey the equations 

y A ' =M A (A,C)y A , z A ' = -z A M A (A, C), (4.2) 

where 

\ 



' 

1 

The eigenvalues of A4^ a (X) = A4 A (X, ±oo) are given by /i Q (A) + Hp(X). The eigen- 
value A* A (A) = /Ji(A) + ^3 (A) with the minimal real part is simple in the right 
complex half-plane Re A > 0. Therefore, the right r A (A) and left 1 A (A) eigenvectors 
of the matrix M. ^(A) associated with this eigenvalue are analytic everywhere in the 
right complex A-half-planc. This, in its turn, implies the analyticity of the vector- 
functions y A (A, C) and z A (A, £) everywhere in the right complex A-half- plane. The 
components of the vectors r A (A) and 1 A (A) are given via the components of the 
vectors ri^A) and li, 3 (A) by the expressions similar to (4.1) everywhere where the 
last vectors are well defined, i. e. in some neighbourhood of the origin of the complex 
A-plane. The vector-solutions y A and z A satisfy 

lim c-^ A ( A V A (A,C) =r A (A), 

lim c^ A V(A,C) =1 A (A). 

The analytic in the complex half-plane Re A > Evans function -D(A) is defined 
as follows: 

z 1 (A,C)-yi(A,C) z 1 (A,C)-y 3 (A,0 \ u ^ 
z 3 (A,C)-yi(A,C) z 3 (A,C)-y 3 (A,C) )' ( } 

The last equality in (4.3) has the sense where the determinant in the right hand side 
is well defined, in particular, in the neighborhood of the origin, where, according 
to Theorem 3.2, the vector-functions yfe(A, £), Zfe(A, (), k — 1,3 are analytic. From 
(4.2) it follows that D(X) is independent on (. It also follows that the function D(X) 
is analytic in the neighborhood of the origin of the complex A-plane. 

Theorem 4.1. ([2], [3]). For Rc A > the function D(X) has zeroes if and only if 
there is a solution of (1.4), which decays exponentially as ( — > ±oo. 



M A (A,C) = 








1 





















1 


1 







2Ac 


-2p ' 


1 - c 2 - p° 










1 


















1 


A 2 


+ P°" 








1-c 2 


-p° 










A 2 + p°" 





-2Ac + 


2p ' 






D(X) = z A • / 



det 
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The asymptotic behaviour of the Evans function D(X) as A — > oo is governed 
by the following 

Theorem 4.2. D(X) -» 1 A (A) • r A (A) as |A| -» oo. 

To prove Theorem 4.2 it is sufficient to verify the conditions of proposition 1.17 
in Pego and Weinstein [20] (we note that in [20] the vectors corresponding to r A (A) 
and 1 A (A) are normalized so that 1 A (A) • r A (A) = 1, and their Evans function has 
asymptotic value unity). It has to be checked that 

• M^ D (X) is diagonalizable for large A; 

• = W(X)H(C)V(X) -> as |A| -> oo, where V(A) is a matrix of right 
eigenvectors, W(A) is its inverse and TZ(() = M A (X, () — M^o(X). 

The validity of the first condition is given by the 

Lemma 4.3. ([2]). For large \X\ the matrix M /K {X) has eigenvalues 

±V2A(1 + OdAp 1 )), iiV^l + OdAp 1 )), ±ic(l + 0(|A|- 2 )). (4.4) 

The second condition is verified by direct computation using Mathematica 4 . 0. 
Consider the eigenvectors v^, i = 1,2,3 associated with the eigenvalues (4.4) with 
first component normalized to unity. One has = v^(l + 0(|A| -1 )), where 

„± / /-— l-c 2 ±cV2A + 2A -1 + c 2 =fcV2A + 2A 
Vl = ^1;±V2A; ; ; 

(l-c 2 ±cV2A-2A)VA 2 cVA(-l + c 2 T\/2Ac + 2A)\ T 

T 7=2 ;A T 7=2 ) ' 

A . . r— 1 - c 2 =F ic\/2X — 2X -l + c 2 ±ic\/2A-2A 
v 2 = I 1; ±iV2A; ; ; 

i(-l + c 2 ±ic\/2A-2A)\/A 2 icVA(-l + c 2 ±i\/2Ac-2A)\ T 

± vf ;A± T2 ) ' 

«± (, , • 1 2 -x -l±2iA -c(±i + 2A) Wl . 
vf = \ l; ±ic; 2~ lA; 2 ; "^71 — ; ~ X{ + } 

Then V = ^(l + OdAI- 1 )), where V is a column matrix combined from the vectors 
v^. The determinant A of V is given by the expression A = A(l + 0(|A| -1 )), where 

A = -16A 2 (c 4 -4A 2 ) 2 

is the determinant of V. Hence, both V and V are invertible for large |A| and T = 
T(l + OdAI" 1 )), where f = WKV and VV = V" 1 . The direct computation shows 
that T = A _1 A/", where the components of the matrix Af(X, () arc polynomials in 
y/X and the maximal degree of these polynomials is 11. Therefore, F{X, Q — ► as 
|A| — > oo and the proposition 1.17 of Pego and Weinstein [20] may be used in our 
case. 
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5. Computation of the leading coefficient of the Taylor series 

of D(X) 

Let us denote Y k the first, and Z k the last component of the vectors and z k , 
obeying at the same time (1.4), (2.4), respectively. From (3.4) and (3.2) one has 



Ks(A,C) = l- 
as — > co, and 

2 = C »T^*< 



1 + A A + 



c 2 C 2 1 + 2c 2 ' 

+ — — C 



1 - c 2 (1 - c) 



2(1 -c 2 ) 2 2(l-c 2 ) 5 /2 
!>A 2 + 0(A 3 ); 



A 2 )+0(A 3 ), 
(5.1) 



1 -c 2 



c 2 C 2 



l + 2c 2 



2(1 -c 2 ) 2 2(l-c 2 ) 5 /2 



A 2 + 0(A 3 ), 



C 



1 c 2 



Z, i A.u^-+ T ^A + -^- F A 2 + 0(A 3 ), 



(5.2) 



as ( — > — oo. 

According to Theorem 3.2 we look for solutions of the equations (1.4), (2.4) in 
the neighbourhood of A = in the form of expansions 

Y k (X) = Y k0 + XY kl + X -\ 2 Y k2 + 0(A 3 ), 
Z k {\) = Z k0 + XZ kl + l -\ 2 Z k2 + 0(A 3 ). 
The coefficients of the above expansions obey the equations 

CY ko =0, (5.3) 



C* Z, 



kO — 0, 



(5.4) 



CY kl =2cF fc ' , 



(5.5) 



C* Z kl = -2cZ' k 



fcO> 



CY k2 = -2Y k0 +4cYl 1 , 



(5.6) 
(5.7) 



C*Z k2 = -2Z k0 -4cZ' kl , (5.8) 

where the differential operators C and C* are defined as follows: 

LY = Y"" - (1 - c 2 )Y" + (p°Y)", C*Z = Z"" - (1 - c 2 )Z" + P °Z". 

The equations (5.3-5.8) have to be solved taking into account the conditions (5.1), 
(5.2). 
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For A = the complete bases of solutions of the equations (1.4) and (2.4), which 
coincide with the equations (5.3) and (5.4), are determined explicitly. Putting A = 
in (1.4) and integrating it twice we get 

d?w 

C a w=—^-{\- c 2 )w +p°w = a + bC (5.9) 

The solving of the equation (5.3) is equivalent to solving of the second order equa- 
tion (5.9) with arbitrary constants a and b. When a = b = the particular solutions 
are 

W\ = , W2 — Wl J w^ 2 d(. 
For some a ^ 0, b = the solution of (5.9) is given by 

w 3 = -wi J w 2 d( + w 2 J wi d(. 

For a = and some b ^ the solution of (5.9) is 

w A = -wi j (w 2 d( + w 2 J (wi d(. 

The solutions wi, w 2 , u> 3 , w 4 constitute the basis in the space of solutions of 
the equation (5.3). By the appropriate normalization of wi and W3 one gets Y10 
and F30 correspondingly, such that according to (5.1) Yio — » and Y30 — > 1 as 

( — > 00: 

e sC - e-< _ 1 - 6c 2 < + e 4 < 

10 ~ (cX+e-<) 2 ' 30 ~ (l+c 2 <) 2 ' 

where s = \/i — c 2 hereafter. 

Normalizing W2 one gets Y 2 o, having asymptotic Y 2 o — > e sl » as £ — > 00. Linear 
combination of W4, Y03 and T20 gives I40 asymptotic to ( as £ — *• 00. 

Putting A = in (2.4) we get the equation (5.4). Following [2] we denote w*" = 
4>. Then <f> = wi and <p — w 2 are two solutions of (5.9) with a = b = 0. Therefore, 

<= J (yj T 2 d c) <%> w* 2 = y (y ^ 

satisfy the equation (5.4). Normalizing to get the asymptotic (5.2) we obtain 

Zio = arctane* 1 ', 

Zio — > e s< * as £ — > 00. The solution Z 2 o asymptotic to c - ^ as ( — > 00 is got by 
normalizing ^2 • The solutions Z 30 = 1 and Z40 = C complete the basis in the space 
of solutions of the equation (5.4). 

The equations (5.5), (5.7) have the form CY = W. Integrating the last equation 
twice one gets 



/(/ 



C Q Y = W, W= / / Wd()d(. (5.10) 
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The particular solution of (5.10) is given by 

Y = -Y w J WY 2Q d( + y 20 J WY W dC (5.11) 

The symbolic language Mathematica 4.0 was employed for construction of the 
solutions of the form (5.11) in the explicit form. To get the asymptotics (5.1) we 
modify the solutions of (5.5), (5.7) in the form (5.11) by adding a linear combination 
of the solutions Y a0 of the homogeneous equation (5.3). For example, 

To get Yii we add 7F10 to the solution of (5.5) in the form (5.11), where the constant 
7 is determined by (5.1). 

The equations (5.6), (5.8) have the form C*Z = W. Putting <j) — Z we get 
the equation of the form (5.10) for the function <fi. Using the formula (5.11) to 
obtain the particular solution, we come back to Z integrating <p twice. Again, to 
get the asymptotics (5.1) wc modify the solutions of (5.6), (5.8) by adding a linear 
combination of the solutions Z a Q of the homogeneous equation. For example, 

carctan(e s< ») c£ arctan(e sl >) icPolyLog [2, — ie s ^} ic PolyLog [2, ie sl »] 
Zn = ^ ■ 7. 1 



^31 — , Z; 



c' ~ 32 ~ s 4 (l-c) 

rCe sC , Z 32 ^-^^ 7 C^-^ 



2s 3 2s 3 
s 2 (l + c)C 2 -4(l + c) ln(l + e 2 <) 



where PolyLog [2, y] = J° t _1 ln(l — t) dt, so that 

i PolyLog [2, -ie sC ] - i PolyLog [2, ie sC ] -> 2e sC , C -» -00, 

i PolyLog [2, -ie sC ] - i PolyLog [2, ie sC ] tt,< + 2e" sC , (^00. 

To get Zn and Z 32 we modify the solution of (5.6), (5.8), obtained with the help 
of (5.11) and further integration, by adding 71Z10 and 72^40, respectively, for the 
constants 71,2 to be determined from (5.2). 

The expressions for the other solutions of (5.3)-(5.8) are rather long and they 
are given in Appendix A. 

The coefficients at the powers of A for the expansion of the vector functions 
yi,3(^)C); z i,3(XC) are uniquely determined by the corresponding functions Y^, 
Z%k, i = 1) 3, k = 0, 1, 2 (see the expressions at the bottom of (2.3) and (2.5)). 

To get the coefficients of the expansion series in a zero neighbourhood for D(X) 
the formula (4.3) was used. We compute 

f\(\\ — A f f fl H + 6 H A + d " A2 + ( A3 ) fl 12 + b ^ X + + C(A 3 ) 

[ ' { 02i + &21A + d 21 X 2 + 0(A 3 ) a 22 + b 22 \ + d 22 \ 2 + 0(A 3 ) 



where 



a m „ = 0, m,n= 1,2, 6 n = b 2 i = d 2 i = 0, 6 12 = — — — -, 6 22 = 2, 

2(1 + c) 

l-2c 2 2 + 3c 4(1 + c) 

dll = ^ 3- ' rfl2 = ^2(iT^' d22 = "(T~^- 
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Since a mn = the first nonzero coefficient of the expansion of D(X) is that one at 
the third power of A and, moreover, this coefficient is completely determined by the 
coefficients at the first and second powers of A of the elements of the determinant. 
Finally, we have 

in some neighbourhood Uq of the origin in the A-plane. 

According to Theorem 4.2 the normalized Evans function D(X), D(X) — > 1 as 
|A| — > oo, is given by 

° {X) ~ 1A(A)tA(A)- 
In a zero neighborhood one has 

^-"UK! !M)--*-^ +0 <*>- 

Hence, we get (1.5) and Theorem 1.2 is proved. 

The computation of the normalized Evans function D(X) can be also done by 
means of the integral formula (Pego and Weinstein [20], Theorem 1.11), which in 
our case has the form 



oo 



D(X) 



J z-(A,C)-Al^(A,C)-y A (A,C)ciC, (5.12) 



where dot denotes differentiation with respect to A. The straightforward computa- 
tion gives 

D{X) = -2 l -^X + 0{X% 

that is in full agreement with (1.5). We list the expresion for the coefficient at A 2 
of the integrand in (5.12) in Appendix A. 



6. Conclusion and Discussion 

In this paper we established that for c 2 < 1/2 the leading coefficient of the Taylor 
series of the normalized Evans function D(X) in a neighbourhood of zero is negative 
and thus, D(X) has to vanish somewhere on the positive real axis. It follows then 
from Theorem 4.1, that there exist the unstable eigenvalue and associated with 
it unstable eigenfunction given by (1.3). This in its turn, implies the exponential 
instability of the loop solitary wave in the prescribed velocity range. 

In all cases known to us, when the equation under analysis has a hamiltonian 
structure and an invariant momentum associated with translational invariance of 
the problem, there is the link between the convexity of the momentum and the sign 
of the leading Taylor coefficient (Pego and Weinstein [20], Alexander and Sachs 
[2], Bridges and Derks [7]). In these cases the existence of the unstable eigenfunc- 
tion is caused by the translational symmetry. In the case under consideration the 
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momentum associated with translational invariance T(w)w(£) = expfwc^) w(£) = 
w(C + oj), u) G R, is Q = Q(w) (see Introduction). One has 

jgW _ 8 c Q 
dc s 3 

i.e. the momentum is concave for all c e [0,1). Yet, the lead coefficient in the 
expansion of D(\) changes its sign at c = l/\/2 and, therefore, the above rule about 
the link is violated. The reason of this violation is that the other symmetry than the 
translational one is responsible for the existence of the unstable eigenfunction in 
our case, namely the rotational symmetry around the x\ axis: Q(ip)w = cxp(Aip) w, 
cp E Ei 1 , where A is the block diagonal 6x6 matrix with the blocks 




The question about stability of the loop solitary wave for c e [1/V2, 1) is still 
open, however it looks plausible that the impulse of the moving loop for this range 
of speeds is large enough to stabilize it and the loop manifests at least marginal 
stability. This conjecture is based on the analogy with a different problem of stability 
of solitary impulses in a composite medium. The linearized problem about the 
solitary impulse in the composite is exactly the same as the linear problem (2.2) of 
the paper, and it is known from the rigorous nonlinear stability analysis (and also 
from numerical evidence) that solitary impulses in question are orbitally stable. It 
may mean that there is no unstable eigenfunction for (2.2) for c > l/\/2. 

Let us explain this analogy in more detail. The system of governing equations, 
describing propagation of quasi-transverse elastic waves in some model of composite 
with nonlinear anisotropic elastic matrix has the form (see e.g. Il'ichev [14]) 

dui dv t dvi d ( d@ \ d 3 u t . 



dt* dx ' 9i* dx \ dui J dx 3 

Here po is the average density, Uj, Vi - gradients of displacements and particle 
velocities in the wave front, respectively, <& is the elastic potential, given by the 
expression 

* = \f{u\ + u\) + \g{u\ - u\) - \k(u\ + ulf. 

The constants / > 0, g > 0, k and m > characterize elastic modules, anisotropy, 
nonlinearity and dispersion, respectively. There is one-parametric anisotropic family 
of solitary wave solutions of (6.1) for k > 

u? = ± y /2 P0 K- 1 (iJL 1 - V 2 )scch^ pam-^m - V 2 )(x - VU), v% = 0, (6.2) 

where pi = (f — g)/po and V is the dimensional speed [14]. The solitary wave 
orbit is stable for pi/2 < V 2 < pi ([14], Theorems 4, 5). The results of numerical 
calculations reported in Bakholdin et al [5] confirm this theoretical result, and also 
give the evidence that solitary wave family (6.2) is unstable for < V 2 < pi/2. 
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The equations (6.1), being linearized about the solitary wave (6.2) for small 
perturbations 8u\ =u\ — u\, arc transformed into one linear equation 

d 2 5ui d 2 Sui k , r ) \id' 2 5ui m d 4 6ui , 

-Mf = ^^- 3 7o {ui) " 7o~w- (6 - 3) 

Next, define the dimensionless speed c e [0, 1) and variables £, t, r by 



c=-^=, x = J-^—S,, t* = S-^jt, 5ll\ = \l———T. 

Now, in the dimensionless variables, we have that solitary pulses in the composite 
are nonlinearly orbitally stable for c 2 > 1/2 (Theorems 4,5 of [14]). In these variables 
the linearized problem (6.3) has exactly the same form as the basic linear problem 
for the loop stability (2.2) with p° given by (1.2). This result strongly prompts that 
the linearized problem (2.2) can have no unstable eigenfunctions for c 2 > 1/2, i.e. 
that our loop solitary wave is at least marginally stable inside this range of speeds. 

The present work is fulfilled under the financial support of the Russian Foundation for 
Basic Research grant No. 02-01-00729 and the President Program of Support of Leading 
Scientific Schools grant No. 1697.2003.1. 

Appendix A. Mathematica expressions for Y 3 i, Yi 2 , Y 32 , Z i2 

and for the integrand in (5.12) 

In the expressions below the notation x instead of Q is used. 

Y 31 = -(e- s x (32 c e s x + 32 c 2 e s x - 96 c e 3 s x - 96 c 2 e 3 s x + 2 e s x s 2 - 
2e 5sx s 2 +n s 2 - e 2 s x ir s 2 - 17 e 4 s x tt s 2 + e 6 s x it s 2 +8 c e s x s x+ 
8c 2 e SI si-48ce 3sl s^48c 2 e 3sl si + 8ce 5 "sa; + 
g c 2 e 5 s x s x + 4 e s x s 3 x - 24 e 3 s x s 3 x + 4 e 5 s x s 3 x - 
12 e 2 s x tt s 3 x + 12 e 4 s x tt s 3 x - 2 s 2 arctan[e s x ] + 
18 e 2 s x s 2 arctan[e s x ] + 18 e 4 s x s 2 arctan[e s x ] - 
2e 6sx s 2 arctan[e s x ] + 12 i e 2 s x s 2 PolyLog[2, -i e s x ] - 
12 i e 4 s x s 2 PolyLog[2, -i e s x ] - 12 % e 2 s x s 2 PolyLog[2,i e s x ] + 

12 % e 4 s x s 2 PolyLog[2, i e s x ])) / 

(4 (1 + c) (l + e 2sx ) 2 s 3 ); 

Y 12 = ^{e- s x (c 2 -17 C 2 e 2 "- S 2 + 13e 2 " S 2 - 

2 (1 + e 2 s x f s eV V 

4e 4sx s 2 + 2 e sx it s 2 - 12 e 3sx tt s 2 + 2 e 5sx tt s 2 - 

14 c 2 e 2 s 1 s a; + 6 c 2 e 4 s x s x - 2 e 2 s x s 3 x + 2 e 4 s x s 3 x - 

2 c 2 e 2 s x s 2 x 2 + 2 c 2 e 4 s x s 2 x 2 - 4 e s x s 2 arctan[e s x ] + 

24e 3sx s 2 arctan[e s x ] - 4 e 5 s x s 2 arctan[e s x ] - 

2 s x 2 i „ f-i i —2 s x~\ | o 4 s x 2 f-i , —2 s x 



e s 



\og[l + e- 2sx ]+8e 4sx s 2 log [l + e~ 2 s x ] )) ; 
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Y 32 = , (e~ s x ( - 128 c 2 e s x - 128 c 3 e s x + 

2 (1 + c) (1 + e 2 **) 2 S 6 V v 

256 c 2 e 3sl + 256 c 3 e 3 " - 2 e s 1 s 2 + 6 c e" s 2 - 52 e 3 s 1 s 2 - 

124 c e 3 s 11 s 2 + 6 e 5 s x s 2 + 14 c e 5 s x s 2 - 2 tt s 2 + c tt s 2 + 

2 e 2 s x tt s 2 - 10 c e 2 s x tt s 2 + 34 e 4 s x tt s 2 + 62 c e 4 s x tt s 2 - 

2 e* 3 s x tt s 2 — 5 c e* 3 s x tt s 2 — 64 c 2 c s x s x — 64 c 3 c s x s x -\- 

192c 2 e 3sx sx + 192 c 3 e 3sx sx- 18 e sx s 3 x- 50c e sx s 3 x + 

96 e 3 s x s 3 x + 192 c e 3 s x s 3 x-U e 5 s x s 3 x - U ce 5 s x s 3 x + 

c tt s 3 x + 24 e 2 s x tt s 3 x + 11 c e 2 s x tt s 3 x - 24 e 4 s x tt s 3 x - 

53ce 4sx tt s 3 x + ce 6sx tt s 3 x-8c 2 e sx s 2 x 2 -8c 3 e sx s 2 x 2 + 

48 c 2 e 3sx s 2 x 2 + 48 c 3 e 3 s x s 2 x 2 - 8 c 2 e 5 s x s 2 x 2 - 8 c 3 e 5 s x s 2 x 2 - 

2 e s x s 4 x 2 - 6 c e s x s 4 x 2 + 12 e 3 s x s 4 x 2 + 36 c e 3 s x s 4 x 2 - 

2e 5 s x s 4 x 2 -6 ce 5 8 x s 4 x 2 -12ce 2 s x tt s 4 x 2 + 12ce 4 s x tt s 4 x 2 + 

4 s 2 arctan[e s x ] - 2 c s 2 arctan[e s x ] - 100 e 2 s x s 2 arctan[e s x ] + 

22ce 2sx s 2 arctan[e s x ] + 28 e 4 s x s 2 arctan[e s x ] - 

94ce 4sx s 2 arctan[e s I ]|4e 6sl s 2 arctan[e s x ] + 

10 c e 6 s x s 2 arctan[e s x ] + 2 s 3 x arctan[e s x ] + 

2 c s 3 x arctan[e s x ] - 18 e 2 s x s 3 x arctan[e s x ] - 

18 c e 2 8 x s 3 x arctan[e s x ] - 18 e 4 s x s 3 x arctan[e s x ] - 

18ce 4 "^a: arctan[e s x ] + 2 e 6 8 x s 3 x arctan[e s x ] + 

2ce 6sx s 3 x arctan[e s x ] - 16 c e s x s 2 log [l + e~ 2 s x ] + 

96 c e 3 s x s 2 log [1 + e- 2sx ]-16ce 5sx s 2 log [l + e" 2 s x ] + 

8e sx s 2 log [l + e 2 "]+8ce" S 2 log [l + e 2 s x ] - 



3 s x s 2 log [1 + e 2 s x ] - 48 c e 3 s x s 2 log [l + e 2 s x ] + 
5 s x s 2 log [1 + e 2 s x ] + 8 c e 5 s x s 2 log [l + e 2 s x ] - 



8e 5sx s 2 log [1 + e 2 s x ] + 8 c e 5 s x s 2 log [l 

1 s 2 PolyLog[2, -ie sx ]-2ics 2 PolyLog[2, -i e s x ] - 15 i e s 

PolyLog[2, -ie sx ]+%ice 2sx s 2 PolyLog[2, -i e s x ] + 
33 i e 4 s x s 2 PolyLog[2, -i e s I ]+54ice 4 " s 2 PolyLog[2, -i e s x ] - 
ie 6sx s 2 PolyLog[2, -i e s x ] ~ 2 i c e 6 s x s 2 PolyLog[2, -i e s x ] - 
12ie 2sx s 3 x PolyLog[2, -i e s x ] - 12 i c e 2 s x s 3 x PolyLog[2, -i e s x ] + 

12 i e 4 s x s 3 x PolyLog[2, -i e s x ] + 
12ice 4sx s 3 x PolyLog[2, -i e s x ] + i s 2 PolyLog[2, i e s x ] + 
2ics 2 PolyLog[2, i e s x ] + 15 i e 2 s x s 2 PolyLog[2, ie sx }- 
6ice 2sx s 2 PolyLog[2, i e s x ] - 33 i e 4 s x s 2 PolyLog[2, i e s x ] - 
54 i c e 4 s x s 2 PolyLog[2, i e s x ] + i e 6 s x s 2 PolyLog[2, ie sx ] + 

2 i c e 6 s x s 2 PolyLog[2, i e s x ] + 12 i e 2 s x s 3 x PolyLog[2, ie sx ] + 

12ice 2 sx s 3 x PolyLog[2, ie sx ]- 
12ie 4sx s 3 x PolyLog[2, i e s x ] - 12 i c e 4 s x s 3 
x PolyLog[2,i e sx ] + 24ie 2sx s 2 PolyLog[3, -i e s x ] + 
48 i c e 2 s x s 2 PolyLog[3, -i e sx ]-2Aie 4sx s 2 
PolyLog[3, -ie sx ]-48ice 4sx s 2 PolyLog[3, -i e s x ] - 

24 i e 2 s x s 2 PolyLog[3, i e s x ] - 
48ice 2sx s 2 PolyLog[3, i e s x ] + 24 i e 4 s x s 2 PolyLog[3, i e s x ] + 

48 i c e 4 s x s 2 PolyLog[3,i e s x })); 
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Z 12 = ( - c 2 e" x + e s x s 2 + 7 c 2 arctan[e s x ] + 

13 s 2 arctan[e s x ]+2 c 2 s x arctan[e s x ] + 2 s 3 x arctan[e s x ] - 
2 c 2 s 2 x 2 arctan[e s x ] - 16 s 2 arctan[e s x ] log[2] - 
i c 2 PolyLog[2, -i e s x ] - 5 i s 2 PolyLog[2, -i e s x ] + 
2 i c 2 s x PolyLog[2, -i e s x ] + i c 2 PolyLog[2, ie sx ] + 
5 i s 2 PolyLog[2, i e s x ] - 2 i c 2 s x PolyLog[2, ie sx ] + 

8 i s 2 PolyLog [2, 1 - 1 i e s *] - 8 i s 2 PolyLog [2, i + 1 i e s x ] - 

2 i c 2 PolyLog[3, -ie sx ]-2i s 2 PolyLog[3, -i e s x ] + 

2 i c 2 PolyLog[3, ie sx ]+2i s 2 PolyLog[3, i e s x ]) ; 



Integrand = (e s x ( - 16 c 3 e s x + 16 c 5 e s x + 16 c 3 e 3 3 x - 
16 c 5 e 3 s x + 16 c e s x s 2 - 16 c 3 e s x s 2 - 16 c e 3 s x s 2 + 
16 c 3 e 3 s x s 2 + 2 c TT s 2 - 2 c 2 tt s 2 - 12 c e 2 s x tt s 2 + 
12 c 3 e 2 s x TT s 2 + 2 c e 4 s 11 tt s 2 + 2 c 2 e 4 s x tt s 2 - 4 c 3 e 4 s x tt s 2 + 
tt s 4 - c tt s 4 - e 4 s x tt s 4 + c e 4 s 11 tt s 4 - c 2 tt s 3 a; + c 3 TT s 3 X + 
6 c 2 e 2 s 35 tt s 3 x - 6 c 3 e 2 s x tt s 3 x - c 2 e 4 s x it s 3 x + c 3 e 4 s 35 tt s 3 x - 
c it s 5 x + 6 c e 2 s x it s 5 x - c e 4 s x ir s 5 x - 8 c s 2 arctan[e s x ] - 
8 c 2 s 2 arctan[e s x ] - 8 c 3 s 2 arctan[e s x ] + 48 c e 2 s x s 2 arctan[e s x ] + 
96 c 2 e 2sx s 2 arctan[e s x ] + 48 c 3 e 2 3 x s 2 arctan[e s x ] - 
8ce isx s 2 arctan[e s x ] - 24 c 2 e 4 s x s 2 arctan[e s x ] - 
8 c 3 e 4sx s 2 arctan[e s I ]-4s 4 arctan[e s x ] + 
4e 4 "s 4 arctan[e s x ] + 2 i c 2 s 2 PolyLog[2, -i e s x ] - 
12 i c 2 e 2 s x s 2 PolyLog[2, -i e s x ] + 2 i c 2 e 4 s 1 s 2 PolyLog[2, -i e s x ] - 
2ic 2 s 2 PolyLog[2, i e s x ] + 12 i c 2 e 2 s x s 2 PolyLog[2, ie sx ]- 

2 i c 2 e 4 s x s 2 PolyLog[2, i e s x })) / 

((-1 + c) (1 + c) (1 + e 2 ^) 3 s 2 ); 



The rest of solutions of (5.3)-(5.8) with their asymptotics are given in sec. 5. The 
listed functions have the asymptotics 



c 2 x 2 1 + 2c 2 



(1-C 2 ) 2 (l-C 2 ) 5 /2- 

x Y ^ x 2 
1-c' 32 (1-c) 2 



5^31^-^ :, Y 3 2^ — T2, as i ^ oo, 



and 



c 2 a; 2 1 + 2c 2 . 

; 1 - " " I (1 _ C 2)2 - (1 _ C 2 )5 /2 X ) e ' aS ' " ^ 
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The function P(x) — PolyLog[3, ie x ] above is denned via 



dP{x) 



= PolyLog[2,ie a; ], 



dx 



and has the asymptotic behaviour 



ImP(ar) -» -e x , 



as x 



— oo; ImP(a;) — > — —a; 2 



as x 



oo. 
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